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Abstract 

We study two aspects of the CFL phase of QCD in the NJL approximation. The 
first one is the issue of the dependence on u of the ultraviolet cutoff in the gap 
equation, which is solved allowing a running coupling constant. The second one 
is the dependence of the gap on the strange quark mass; using the high density 
effective theory we perform an expansion in the parameter (m s /fj,) 2 after checking 
that its numerical validity is very good already at first order. 



1 Introduction 



The existence of color superconductivity at very large densities and low temperature is 
an established consequence of QCD (for general reviews see PQ, and j2j). Since at lower 
1 On leave from the Dipartimento di Fisica, Universita' di Firenze, F50019 Firenze, Italia 
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densities one cannot employ perturbative QCD, a popular approach in this regime is 
the use of a Nambu-Jona Lasinio (NJL) model. The mathematical procedure consists in 
solving the mean field gap equation and selecting the solution which minimizes the free 
energy. At sufficiently high densities, for three massless flavors, the condensation pattern 
leads to a conserved diagonal subgroup of color plus flavor (color-flavor locked phase, 
briefly CFL phase [3 Hj ) ; this phase continues to exist when the strange quark mass is 
non vanishing, but not too large jSlE]- An interesting property of the CFL phase, proved 
by Rajagopal and Wilczek jjj, is its electrical neutrality, which implies that the densities 
of u, d and s quarks are equal. A stable bulk requires not only e.m. neutrality, but also 
color neutrality. It should also be a color singlet, but according to ref. [S] in a color 
neutral macroscopic system imposing this property does not essentially change the free 
energy. Alford and Rajagopal 9j have shown that in neutral CFL (with finite strange 
quark mass) quarks pair with a unique common Fermi momentum. The neutrality result 
is basic to our calculations below. It allows the use of a well defined approximation of 
QCD at high density (high density effective theory HDET), see [TU1 HU H2] and, for a 
review, |13|. 

In this letter we address two aspects of CFL for QCD modeled by a NJL four fermion 
interaction. Even though NJL is only a model, it offers simple expressions that can 
be helpful in clarifying physical issues; therefore a better understanding of its dynamics 
is significant. The two aspects concern the role of the ultraviolet cutoff in the NJL 
interaction and the relevance of the effects due to the strange quark mass. As for the first 
point, the cutoff is usually fixed once for all and considered among the parameters of the 
model. However, when working at varying large densities the choice of the appropriate 
cutoff is rather subtle. It is suggestive to consider, for instance, a solid where there is 
naturally a maximum frequency, the Debye frequency. One expects that when particles 
become closer the ultraviolet cutoff extends to larger momenta. The appropriate physical 
simulation of the real situation would then require a cutoff increasing with the chemical 
potential, rather than a unique fixed cutoff. Our analysis suggests that, in order to get 
sensible results from the NJL model, the ultraviolet cutoff should similarly increase with 
density. The results we obtain for the physical quantities, below in this paper, show 
unequivocally that this is indeed the case. In the text we will explain in more detail the 
essential difficulties (such as a decreasing gap for larger densities) in which the theory 
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would run into if taken with a fixed cutoff. We shall propose and apply a convenient 
procedure to solve the problem in terms of a redefinition of the NJL coupling constant 
such as to make it cutoff dependent. In Section 2 we discuss this issue in the CFL model 
with massless quarks and we find the optimal choice for the dependence of the ultraviolet 
cutoff on the quark chemical potential. 

The second aspect we want to discuss is the role played by the non vanishing strange 
quark mass. We address it in Section 3, where we provide semi- analytical results for 
the dependence of the various gap parameters on m s . CFL with a massive strange quark 
represents a more realistic case in which our previous discussion can be applied; we include 
both the triplet and the sextet gaps and perform a perturbative expansion in m^//i 2 , 
obtaining simple expressions for the first non trivial term in the expansion. We compare 
this expansion with the numerical results from the complete gap equations and we find 
that indeed the first perturbative term adequately describes the full dependence. We 
implement from the very beginning the electrical neutrality for the CFL phase and, as 
already observed, this allows the use of the HDET already applied for the 2SC case (Ej; 
for completeness also the simpler massless case is treated by the same formalism. We 
notice that the gap equation with finite strange quark mass has already been discussed 
in 0|T3j. Besides the use of HDET, our main contribution in this context is to provide 
semi-analytical results for the mass dependence of the CFL gaps. We conclude the paper 
with an Appendix where we list some results and integrals related to the gap equation. 



2 NJL running coupling constant 



When the NJL interaction is used for modelling QCD at vanishing temperature and 
density, one can fix the UV cutoff A such as to get realistic quark constituent masses. 
Typically the cutoff is chosen between 600 and 1000 MeV for masses ranging between 200 
and 400 MeV. In any case A is thought of as fixed once for all. This gives no problems 
at zero density, however leads to difficulties when one tries to simulate QCD at finite 
chemical potential. In fact, at finite density one takes as relevant degrees of freedom all 
the fermions with momenta in a shell around the Fermi surface. The thickness of the 
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shell is measured by a cutoff 5, which is the cutoff for momenta measured from the Fermi 
surface. This cutoff is chosen to be much smaller than the chemical potential p and much 
larger than the gap; S is related to the NJL cutoff A by the relation A = p + 5, because A is 
the greatest momentum allowed by the NJL model. This relation is however problematic 
when one is interested in the behavior of the theory for varying p. The constraint A = p+5 
would force 6 to vanish for increasing p, starting from p < A. In turns this gives rise to 
a vanishing gap. In fact, we recall from the simplest version of the BCS theory that the 
gap has a typical behavior A m 25exp[— 2/{Gp)\, where p is the density of the states at 
the Fermi surface (p oc p 2 ) and G is the NJL coupling constant. Therefore decreasing 
the volume of the shell has the effect of reducing the gap, with a quantitatively different 
reduction from the state density p and the thickness 5. The decreasing of A with p does 
not correspond to the asymptotic (p — > oo) QCD behavior, which is characterized by 
an increase of the gap with p, though with a vanishing ratio A/p 16.. The uncorrect 
behavior of A arises because the model is taken to be valid only for momenta up to A 
which forbids to go to values of p of the order or higher than A. Clearly this constitutes 
an obstacle in physical situations where the typical chemical potential is about 400 or 
500 MeV (e.g. in compact stellar objects) with a 5 of the order 150 or 200 MeV. In fact 
it turns out difficult, if not impossible, to explore higher values of p for any reasonable 
choice of A. 

In this paper we make a proposal to overcome this situation. We start by noticing 
that the phenomenology at zero temperature and density 2 is completely determined by the 
meson decay coupling constant /„- and by the constituent quark mass, or equivalently by 
the chiral condensate. In the NJL model these two quantities are fixed by two equations 
(see for a review JZj), one fixing and the other the chiral gap equation. These equations 
depend on the cutoff A and on the NJL coupling G. Our proposal consists in assuming 
the coupling G to be a function of A in such a way that /„- assumes its experimental value 
and that the chiral gap equation is satisfied for any choice of A. 

To be more explicit, we write the Nambu-Jona Lasinio equations with a three dimen- 

2 For the sake of discussion we consider here the ideal case of massless quarks, the more realistic case 
of a massive strange quark will be considered in the following 
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sional cutoff A [T7]; the equation for the ir leptonic decay constant reads : 



3m 



* 2 



(2tt)3 m 



2tt 2 



A , A 

— arcsmh — 

Vm* 2 + A 2 m* 



where m* is the constituent mass at // = which is determined by the self-consistency 
condition: 

Am* M p 2 dp 

(2) 



m 



4m* ... /" A 



m is the quark current mass which is assumed in this Section to be zero. G(A) is the 
NJL coupling having dimension mass -2 ; it could be understood as the effect of a fictitious 
gluon propagator [3]: 

gfivfiab 

1 D ab = 1 > ( 3 ) 

and A 2 G(A) would take the role of the square of the strong coupling constant. Eq. (JTJ), 
with f n = 93 MeV, implicitly defines the function m* = m*(A) which we use in eq. (0) 
to get the function G = G(A). The result of this analysis is in fig. 1. 
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Fi gure 1: The running NJL coupling constant G (dashed line) and the running constituent mass m* 
(in MeV, solid line) as functions of A (in MeV). A is the ultraviolet cutoff. The vertical axis on the left 
refers to A 2 G(A), while the axis on the right refers to m*. 



Our choice implies that a NJL model is defined at any scale by using the appropriate 
G(A). Whereas in the usual case we have to keep the momenta smaller than the cutoff, 
now, for any given momenta, we can fix the cutoff in such a way that it is much bigger 
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than the momenta. The phenomenology of the chiral world is clearly unaffected by this 
procedure; in fact it turns out that the constituent mass acquires a weak dependence 
on the cutoff, and therefore it can be fixed at the most convenient value. Also the 
quantity G(A)A 2 decreases weakly with the cutoff. In applying these considerations to 
the calculations at finite density, we have only to use the appropriate value of the coupling 
as given by + 5), where now // + S has nothing to do with the value of A chosen to fit 
the chiral world. To give an explicit example we consider the CFL phase with massless 
quarks. There are two independent gaps A and A 9 (A 9 = — 2A if the pairing is only in 
the antitriplet channel) and the gap equations are [3J [T3] : 

A = — Aq arcsmh-j-r — r — 2 A arcsmh-r— 7 , 

6tt 2 \ |A 9 | |A|/ 

VGA . S 
A 9 = arcsmh— . (4) 

07T Z A 



If one uses a fixed value for A = fx + 5, as for instance in ref. jT8], one gets a non monotonic 
behavior of the gap, as it can be seen from fig. |21 (dashed line); a similar behavior was 
found in ^H] (their fig. 1), albeit a different choice of the parameters produces some 
numerical differences. On the other hand the solid line shows an increasing behavior of 
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Figure 2: The CFL gap A for massless quarks, as obtained from eq. 01 versus the quark chemical 
potential for the two cases discussed in the text. Solid line: running NJL coupling G(fi + 5) and cutoff 
S = c/i, with c = 0.35; dashed line: 5 = A — fi, where A = 800 MeV, and G(A) = 13.3 GeV~ 2 . The 
picture shows the different qualitative behavior with fj, of the gap. 

the gap. We see that in this way one reproduces qualitatively the behavior found in 
QCD for asymptotic chemical potential JE]- This result is obtained by the running NJL 
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coupling + 5), with the following choice of the cutoff 5: 

5 = cfi (5) 

with c a fixed constant (c = 0.35 in fig. |2J. The reason for this choice is that, as discussed 
above, when increasing /i, we do not want to reduce the ratio of the number of the relevant 
degrees of freedom to the volume of the Fermi sphere. Requiring the fractional importance 
to be constant is equivalent to require eq. 

In fig. El we plot the gap parameter Ai (the results for Ag are similar) as a function 
of c for three different values of the chemical potential. In general there exists a window 
of values for c: 

c= 0.35 ±0.10, (6) 

where the gap parameters are less dependent on c. This is therefore the range of c we shall 
assume below. It can be also noted that the result for the QCD superconducting model 

A 1 (MeV) 100 
80 
60 
40 
20 
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Figure 3: Gap parameter Ai as a function of c, where 5 = c/x; the quarks are massless. The three 
curves are obtained for three values of the quark chemical potential: /.i = 1.0, 0.7, 0.5 GeV. 

we are considering here has its counterpart in solid state physics, where the analysis of 
19 shows a linear increase of 5 with [i similar to eq. 

These results, valid for the massless confirmed by a complete numerical 

analysis including the effect of a strange mass. This will be discussed in the subsequent 
section. 
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3 Effective lagrangian for gapless quarks 



The massive case is considerably more involved, and the simple set of equations (jlj has to 
be substituted by a system of 5 equations 5j that can be only solved numerically. In JH] 
the CFL phase with massive strange quark was also considered. In comparison with [3] 
the derivation we present here has the advantage of offering semi-analytical results, thanks 
to an expansion in powers of m s j [i. We principally differ from ref. ^B] for the different 
treatment of the cutoff, as discussed in the previous section, and for the inclusion of pairing 
in both the antitriplet and the sextet color channel. The possibility of a semi-analytical 
treatment rests on the HDET approximation. This effective lagrangian approach was 
extended in jT3] to the 2SC phase with massive quarks and here we treat the three flavor 
case. 

In the HDET one introduces effective velocity dependent fields, corresponding to the 
positive energy solutions of the field equations: ip a ,i,n{.%) where a, i are color and flavor 
indices, n = v/\v\, and v is the quark velocity defined by the equation 

p v = [an + r (7) 

with v 1 * = (0, v). The effective fields ip a> in{x) are expressed in terms of Fourier-transformed 
quark fields ip a> i(p) as follows 

1>a,i*{x) =V+J j^y 4 e'^rt-* $ a>i (p) . (8) 

Here V+ is the positive energy projector, defined, together with V~, by the formula 

V ± = 1±((? -; + X - 70) , (9) 

where Xi = rrii/fi and m, is the mass of the quark having flavor i. We now change the 
color-flavor basis introducing new fields ip^ as follows: 

= E^l ( 10 ) 

A=l V 1 



where Xa with A = 1, . . . , 8 are the usual Gell-mann matrices and Ag = y2/3 Ao- We also 
introduce the Nambu-Gor'kov doublet 
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In the xa basis the lagrangian of the quarks, including the quark-gluon interaction and 
the gap therm can be written in momentum space as follows: 



£ = £ + £ 1 + £ A - (12) 
£o is the kinetic term while £\ describes the quark-gluon interaction. They are given by: 



X X Xa 

n A,B=1 

9 

£ i = W X X Xa 

n A,B=1 

We have introduced the symbols 

V AB -^ = Tr[T A T B V]^, 

1 

1 



AB ^ 


-° 





V AB ■£ ) 


Hab 


A 





ii 



Xb 



(13) 



TTil Am 
n AmB A n 



ttH Am _ ~ rp 

n AmB A n — 



V AB • t = Tr 

Tr[T A T m T B V»] A™ , 

T A T m T B V» 



TaTbV 



Am 



(14) 



In these equations T A = Xa/ V%, A™ is the gluon field, and V 1 denotes the matrix 

I v% o \ 
V£ o 

V o o vfj 



I] 



(15) 



with Vf = (1, Vi) for each flavor i; a similar definition holds for V^, with Vf = (1, —Vi). 
In the limit m u = w one has t> „ = t> ^ = 1 , v s — — x 2 s = \J\ — m 2 s / fi 2 s . 

Let us now turn to the gap term £&. We consider CFL condensation in both the 
antisymmetric 3^ and in the symmetric 65 channels. We assume equal masses (actually 
zero) for the up and down quarks and neglect quark- ant iquark chiral condensates, whose 
contribution is expected to negligible in the very large \i limit. Also the contribution 
from the repulsive 65 channel is expected to be small, but we include it because the gap 
equations are consistent only with condensation in both the 65 and the 3 A channels. The 
condensate we consider is therefore 

< ipai C 75 Vte >~ (A - e al3I e l3l + G i3 - (S ai 5^ + 5 aj 6*)) . (16) 
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The first term on the r.h.s accounts for the condensation in the 3a channel and the second 
one describes condensation in 65 channel. As we assume m u = m<f, we put 



A 



A 



A 



ud 



A 
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(17) 



G, 



G 



du 



G u d — Gdd = G\ 



G,, 



Gds = Gi 



Gc.o = Gn 



which reduces the number of independent gap parameters to five. We stress that we 
impose electrical and color neutrality, which, as shown in ref. [7j, is indeed satisfied in 
the color-flavor locked phase of QCD because in this phase the three light quarks number 
densities are equal n u = rid = n s , with no need for electrons, i.e. n e = 0. As a consequence, 
the Fermi momenta of the three quarks are equal: 



PF,u = PF,d = PF,s = Pf 



(19) 



which, in terms of the quark chemical potentials and Fermi velocities can be written as 



Vu\V u \ = Vd\Vd\ = Hs\Vs 



(20) 



It follows that the wave function of the quark-quark condensate has no dependence on the 
Fermi energies and the condensation can be described in the mean field approximation 
by the following lagrangian term containing only the effective fields ip a i,±'- 



Ay e^eij! + dj {5 m 5^ + 5 aj 6* 



(21) 



Using the Nambu-Gor'kov fields one rewrites (}2"Tj) as follows: 

9 / () 



£a = E E x\ 



ft A,B=1 



where 



and 



75 



75 ® A A b 




Xb 



AB \ 



(22) 



(23) 



k AB \ 



( A 12 I 3 

Ah 

|(4A-A 12 ) ^(A-A 12 ) 

V% 1 a a ^ 2 



V 3^4 (A -A 



(24) 



12 , 



•|(A 12 + 2A) J 
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( G x h 

G 2 h 

d-K^-Ga) ^(3G X - 2G 3 - G 2 ) 



(25) 



V ^(3Gi -2G 3 -G 2 ) §(3Gi + 2G 2 + G 3 ) J 
where I n stays for the nxn identity matrix. The nine eigenvalues of the matrix Aab are 
reported in the Appendix. From the lagrangian 



£o + £a = Xa 



V 



AB 



n A,B=1 \ 75 ® A 

the fermionic propagator is obtained 



t 

AB 



75 ® A A b 



»=EE X^liWXB (26) 

ra A,B=1 



AV-fA^V-l- A 2 

1_ 

V AV-fA^V-l-A 1 



AV • iA- 1 



■75 



'75 



AV-IA^V-l-A 1 



AV 



® A \ 



/ 



AB 



The off-diagonal matrix S\ 2 B (1, 2= NG indices) is the anomalous component of the quark 
propagator, which is what we need to write down the gap equations. In matrix form they 
are as follows 



^AB 



i G(/j + 5) p 2 F 

47T 3 



TT[l TTU* 

n AaC n DbB 



+<5 

di\\ 

—5 J — oo 



di Q 5 ab g^S x c 2 D {l) 



(27) 



where the color-flavor symbols H^ aB have been defined in ()14|) and the running NJL 
coupling is computed at \i + S, according to the previous discussion. In eq. (|2"7)l 8 = c/x 
is the cutoff discussed in section |21 The gap equations can be numerically solved; results 
for m s = 250 MeV and c = 0.35 are reported in Fig. 0] 

A semi-analytical solution can be found by performing an expansion in the strange 
quark mass: 

m s m s 

x s = — ss — < 1 . (28) 
Us ii 

Here \i = fi u = fid, fi s = fi + 8fi with S/j, = 0{rn s j \x) 2 . We define: 

Ai = A12 + Gi = A? + t x x\ , A 9 = -2A 12 + 4Gi = A° 9 + e 9 x 2 s , 



A = 



4A° - A° 2 _ 2A^ + A° 2 _ 2A? + A° 

■ — 4iX s , G 2 = t,2% s , "3 = 



6 



6 



6 



-6*2 ,(29) 
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Figure 4: The five gap parameters vs. the chemical potential (i; results obtained by the numerical 
solution of the exact gap equations, for m s — 250 MeV and 5 = 0.35/z. The upper curve refers to 
A12, while the middle one refers to A. The lower curves are the data obtained for the three sextet gap 
parameters |Gi|, | G2 1 and | C?3 1 . Gaps and /i are expressed in MeV. 

where A^Ag are the values for massless quarks and can be obtained from eqns. (@J). 
For m s 7^ the number of gaps increases from 2 to 5, but eqns. ()29j) give immediately 
the solutions for any value of m s (with the proviso x s = m 2 s /{i 2 <C 1) if one knows the 
parameters £ = (ei, eg, £1, £2, £3)- They can be obtained by solving the system of 5 linear 
algebraic equations 

A • £ = f . (30) 

The matrices A and / are reported in the appendix. Numerical results for the parameters 
£ are in Table 1. For a strange mass of 250 MeV we find the rate A9/A1 equal to —2.365 
for n = 500 MeV and -2.33 for fj, = 1000 MeV. 





ei 


eg 


£1 


£2 


£3 


500 


-35.29 


90.48 


33.63 


-2.06 


-6.20 


700 


-44.86 


102.93 


41.16 


-3.10 


-8.72 


1000 


-61.21 


135.08 


55.69 


-4.42 


-12.23 



Table 1: Values of the £ components, obtained by means of the 5x5 system of linear equations. In the 
gap equations we have made the choice 8 = 0.35^. All the gaps and /x's are expressed in MeV. 
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Figure 5: Gap parameter Ai = A12 + G\ as a function of /i. Left panel: results for two different values 
of the strange quark mass, m s = 150, 250 MeV (curves are obtained by means of the perturbative gap 
equations). Right panel: relative variation JAi/Ai = [Ai(m s ) — Ai(0)]/Ai(0). In this plot 6 = 0.35/1. 
Gaps and jt/'s are expressed in MeV. 

From these results we can compute the dependence of the various gaps on fi. For 
Ai it is reported (for two values of the strange quark mass) in figEJ We stress that the 
approximation we consider here is indeed very good in the range of parameters we have 
considered in the present paper, i.e. m s = 150 — 250 MeV and fi = 500 — 1000 MeV. As 
an example, fig. |H1 shows the approximate solution for Ai (solid line), which differs only 
by a few percent from the exact solution (dashed line), except in the lower region of \i. 
These results are obtained for a mass of 250 MeV; for smaller values the effect is even less 
relevant. 

4 Conclusions 

A main point in this work has been to clarify the issue of the dependence on /i of the 
ultraviolet cutoff in the NJL model and its application to CFL phase of QCD. A convenient 
procedure consists in re-defining the NJL coupling such as to make it cutoff dependent. 
The application to the CFL model with m s 7^ leads, in the approximation of the 
high density effective theory, to an expansion in the parameter (m s /fi) 2 , whose numerical 
validity is very good already at first order. Both triplet and sextet gaps are included. The 
gap parameters are obtained as functions of the chemical potential and their behavior is 
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Figure 6: This diagram shows the difference between the exact numerical calculation of the gap Ai for 
m s = 250 MeV (solid line) and the approximate result obtained through the expansion in m s (dashed 
line). Units are MeV. 

consistent with the expected asymptotic limits. 



Appendix 

Eigenvalues of the gap matrix 

The eigenvalues of the gap matrix for CFL with massive strange quark are as follows: 
Ai = . . . = A 5 = A 12 + Gi, A 6 = A 7 = A + G 2 , A 8/9 = \(x-z± J{x + zf + Ay 2 ) 
where x,y,z are given by 

x = i(4A-A 12 ) + G 1 -^(G 2 -G 3 ), (31) 

V = ^(A-A^ + v^f^-^-^V (32) 



z = ^(A 12 + 2A)-^(3G 1 + 2G 2 + G 3 ). (33) 



Matrices £, A, f 
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We write down explicitly the matrices £, A, f defined in the text 



e = (ci e 9 6 6 &) ; 



(34) 



( a [3 

a [3 

a (3 



\ 7 1 



-§K-Si + 2 .99 - 5ig) 
-1- |(75i +459+519) 
-5- f (551 + 259 + 519) 



|&(-5i + 259 - 519) 

§< 



|^(5i - 59 + 2519) 



f (1 



K 5 5i -5ig)) 



-1 + 1(51 + 459 + 519) 
■ - f (751-259-519) 

f (l-|(5i-5i9)) 



|(45x + 259 + 25 19 ) 
; (-5i +59 + 2519) 



26, 
9 



(1 + 1(5! -519)) 



f45! 



166 „ 
-9-91 



2 1 46 „ 

3 + T^l 



(35) 



where we have defined b = G(fi + 5)/j, 2 /2tc 2 , a — 1 — (5 — ^bhg, 7 = y^i! finally 

/ &(-f + f -4(/i + /o-4/i 9 )^ ^ 

fo (-f + f - ^(5/1 - / 9 - 2/ 19 ; 

6 (A (-2/0 + 2/ 9 °) + ^(Uf, -h- 4/ 19 



V 



& (li72(- 7 A - / 9 °) + tf(T/i + /i 9 



(36) 



/ 



The values of the parameters /_,-, J°, /i fc are listed below. 
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A? arcsinh 



AO 



hk = arcsinh 



A" 



- f . 



In this expression we use the following parametric integrals: 



l\x) 



J </-/ - — ' i , ' = 2 7T / ( —arcsinh 



rS l2l V-lV-l + x 2 

drl 

-8 



I s (x 1 ,x 2 ) = j_ 



dH 



j\x) 



8 

2tt 

i 

8 



D(x) 2 

V ■ I V ■ I + xi x 2 
D( Xl )D{x 2 ) 



+ 



X/ 





5 


arcsinh 






x 2 



+ 



xi(xi + x 2 ) ( 
k k — arcsinh 



X2 J-i \ 



(5 




5 






— arcsinh 




)] 


X 2 




Xl 





j\ Xl ,x g ) = J* <Pl 



q 

8 dH mxf z 



% ir arcsinh 



VS 2 + x 2 



l 7T 



,) D( Xl )D(x 9 ) 2{x\ - xl) 
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28 [Jxj + 8 2 -Jx 2 9 + 8 2 



(37) 
(38) 

(39) 



(40) 
(41) 



+2 x 1 arcsinh 



S_ 

Xi 



2 x 9 arcsinh 



6_ 

Xg 



(42) 
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